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ABSTRACT 


A  simple  formula  in  closed  form  is  proposed  for 
the  thermal  accommodation  coefficient  from  which  the 
accommodation  coefficient  of  a  monatomic  gas-solid 
system  may  be  calculated  from  certain  basic  parameters 
of  the  system.  The  formula  is  obtained  from  consideration 
of  certain  aspects  of  Goodman's  lattice  rheory  of 
accommodation  and  from  conclusions  regarding  certain 
properties  of  available  experimental  data.  Adequate 
agreement  is  obtained  of  the  formula  with  both  the  lattice 
theory  and  the  experimental  data;  in  the  authors'  opinion 
all  the  reliable  data  available  are  included.  Certain 
corrections  are  applied  to  the  lattice  theory,  and  the  gas- 
surface  potential  interaction  parameters  relevant  to  the 
lattice  theory  are  revised. 
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INTRODUCTION 


The  thermal  accommodation  coefficient  (hereafter 
abbreviated  to  ac  or  a)  is  a  measure  of  the  average  efficiency 
of  the  energy  exchange  per  encounter  of  a  gas  molecule  with 
the  solid  at  the  gas-solid  interface.  Experimental  methods 
applied  in  the  measurements  of  a  have  been  adequately 
discussed  elsewhere  (see,  for  example,  the  Survey  Article  by 
one  of  us1)*  Here  we  restrict  discussion  to  a  three  dimensional 
(hereafter,  "n-dimensional"  is  abbreviated  to  ”nD" )  monatomic 
Maxwellian  gas  at  temperat  :e  T  which  shares  an  interface 
with  a  solid  maintained  at  temperature  Tg.  If  T'  is  defined 
by 


T'  =  E ' /2k  (1) 

where  E'  is  the  average  energy  per  reflected  gas  atom,  the 
thermal  accommodation  coefficient  is  defined  by 


T-T 1 
T-T_ 


(2) 


The  purpose  of  this  paper  is  to  propose  a  simple  formula  in 
closed  form  for  a  from  which,  it  is  hoped,  the  ac  of  any 
monatomic  gas-solid  system  may  be  calculated  from  certain  basic 
parameters  of  the  system.  This  formula  is  partly  empirical  in 
origin  because,  although  some  of  its  properties  are  obtained 
from  existing  theory,  its  final  form  rests  on  consideration  of 
certain  properties  of  experimental  data. 
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The  first ,  and  perhaps  best  known,  closed-form  formula 
for  a  was  proposed  by  Baule2  ,  who  considered  the  interaction 
in  terms  of  collisions  of  hard  spheres.  The  accommodation 
coefficient  for  the  impact  of  two  hard  spheres,  one  representing 
an  incident  gas  atom  of  mass  M  and  the  other  an  initially 
stationary  solid  atom  of  mass  m,  is 


a 


4y 

(1+y ) 2 


cos2n 


(3) 


where  n  is  the  angle  between  the  velocity  vector  of  the 
incident  gas  atom  and  the  line  joining  the  centers  of  the 
spheres  at  impact  and  y  is  the  mass  ratio  defined  by 


y  =  M/m  .  (4) 

The  Baule  formula  is  obtained  by  averaging  {3)  over  all  n, 

-  i 

putting  cos/n  =  j  : 

(Baule)  a  =  — — —  .  (5) 

(1+y) 2 

Landau3  treated  the  interaction  in  terms  of  classical 
continuum  theory,  using  an  interaction  potential  V(u)  of 
repulsion  only  between  the  gas  atom  and  the  surface, 


V(u)  =  e  exp(-2icu)  (6) 

where  u  is  the  separation  of  the  gas  atom  from  the  surface  and 
e  is  an  energy  of  arbitrary  magnitude.  When  corrected4'5,  the 
theory  results  in  the  formula 


3 


where  0  is  the  Debye  temperature  of  the  solid. 

Gilbey5  used  a  naive  continuum  model  and  obtained 

a  =  1  -  exp(-4E;  cot?)  (8) 

where,  very  approximately 

tan2e  *  w*-l,  (9) 

(the  approximation  (9)  is  not  used  by  Gilbey) . 

None  of  the  above  formulae  yields  extensive  agreement 
with  experimental  data,  although  each  may  be  valid  under 
particular  sets  of  conditions.  There  are  many  formulae  in  the 
literature  which  are  not  in  closed  form;  the  best  known  is, 
perhaps,  the  formula  of  Devonshire4,5.  However,  these  also 
fail  to  yield  extensive  agreement  with  experiment. 

A  sensible  test  of  any  useful  formula  for  acs  requires 
reliable  aata  on  the  temperature  variation  of  acs  of  at  least 
two  known  gases  on  the  same  known  surface.  In  our  opinion, 
the  most  reliable  data  available  are  those  of  Thomas6  and 
Silvernail7  on  the  temperature  dependence  of  the  acs  of  the 
inert  gases  on  W,  and  we  have  placed  considerable  reliance  on 
these  in  this  work. 
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BASIS  OP  THE  FORMULA 

Our  formula  rosts  heavily  on  results  from  the  classical 
lattice  theory  of  acs6"11  (hereafter  referred  to  as  L.T.) 
and  a  3D  theory  of  surface  scattering,12  both  developed  by  one 
of  us,  «3nd  on  conclusions  regarding  certain  properties  of 
available  experimental  data. 

In  L-T.,  the  incident  gas  is  monatomic  and  is  represented 
by  a  ID  Maxwellian  (thermal)  gas  which  impinges  normally  onto 
the  solid  surface,  all  collisions  being  head-on  with  surface 
atoms?  the  solid  surface  is  assumed  to  be  clean  and  is 
represented  by  the  surface  of  a  simple  3D  lattice  model, 
initially  at  0°K.  It  is  assumed  that  a  is  independent  of 
surface  temperature,  consistent  with  data  of  Thomas  and 
Schofield13  in  the  surface  temperature  range  100  -  300°K  a  > ; 
of  Watt  and  Moreton14  in  the  range  800  -  ISOO^K.  Interaction 
potentials  of  the  Morse  type  are  used  between  the  gas  atoms  and 
the  surface  atoms. 

The  results  of  L.T.  are  presented  in  terms  of  a  parameter 
y(t)  called10/11  the  ’’effective  ac  of  a  single  gas  atom”  and 
defined  by  analogy  with  (2)  as  follows.  Initial  and  final 
effective  temperatures,  t  and  tf  respectively,  of  a  gas  atom 
interacting  with  the  surface  are  defined  by 

t  «  MU  2/2k 
o  ' 

tf  ~  MUf2/2k 


and 


(10a) 

(10b) 


c 


where  UQ  and  U^  are  the  initial  and  final  speeds  respectively 
of  the  gas  atom.  Accordingly,  y(t)  for  a  solid  initially 
at  0°K  is  defined  by 

y (t)  =  1  -  tf/t  (11) 

For  a  ID  Maxwellian  gas  at  temperature  T,  the  normalized 
distribution  of  t  is 


*1D 


(t) 


(12) 


and  it  follows  that  the  ac  is  given  by 


(ID  gas) 


a(T) 


t  y (t)  exp  (-  dt 


(13) 


o 

This  is  the  formula  for  a(T)  used  in  L.T. 

A  typical  y (t)  curve  is  shown  in  Fig.  1;  it  is  determined 
essentially  by  four  parameters:  (i)  the  critical  initial 
effective  temperature,  t  ,  for  trapping,  (ii)  the  effective 
temperature,  tm^n,  at  which  y  (t)  has  a  minimum,  (iii)  the  value 
of  Y  (tmn)  and  (iv)  the  high  temperature  limit,  y(«),  of  y(t). 

Y  (t)  =1  for  t  _<  t  ?  y(t)  '  t  /t  for  values  of  t  just  greater 

rTr"~  v«  v 

than  t  ;  for  larger  t,  y(t)  has  a  minimum  and  than  approaches  the 
c 


"hard  spheres  limit",  y{=°)  =  a(®),  as  t  -»■  ».  An  important  result 
of  L.T.  is  that  it  suggests  that  the  shape  of  the  y(t)  curve 
is  independent  of  the  gas-solid  system.  However,  the  values  of 
the  four  parameters  of  the  Y(t)  curve,  t„,  t,.  ,  y(t  .  }  and  y («) , 
depend  on  the  four  parameters  of  the  gas-solid  system  which  enter 


JSX-'SZtf.  • 
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into  the  calculation  of  y(t),  viz.,  M,  m,  0/a  and  D,  where 
a  and  D  are  parameters  of  the  Morse  potential, 

L.T.  considers  only  normal,  head-on,  impact  of  a  gas 
atom  with  a  surface  atom,  and  if  it  is  generalized  to  include 
any  type  of  impact,  a  more  general  y  is  defined.  Accordingly, 
taking  Tg  «=  0,  we  define  y  (t ,  Q ,  4-  ,x ,y)  ,  where  6  is  the  angle 
made  by  the  velocity  vector  of  the  incident  gas  atom  with  the 
inward  normal  to  the  surface,  x  and  y  are  the  coordinates  of 
the  initial  "aiming  point"  on  the  lattice  surface  (we  may  take 
x  =  y  =  0  at  the  center  of  one  of  the  surface  atoms)  and  4>  is 
the  incident  azimuthal  angle  (we  take  <f>  =  0  as  the  x-direction)  . 
For  a  3D  Maxwellian  gas  at  temperature  T,  the  normalized 
distribution  of  t,e,4>,x  and  y  is 

g,n(t,9,<fr,x#y)  =  — - —  sin(2e)  t  exp  [-  (14) 

2  it  AT  2 

where  A  is  the  area  of  the  solid  surface.  It  follows  that  the 
ac  for  a  3D  Maxwellian  gas  is 


(3D  gas)  a(T)  =  — - —  /////  t2  r  (t,  8  ,  <j>  ,x,y)  sin  (26)  exp( -~)dtded$dxdy . 

4  it  AT  3 


(15) 

We  assume  that  the  shape  of  the  r(t)  curve  for  a  3D  gas 
is  also  universal  and  that  in  addition  r(t)  is  a  sum  of  two 
simple  and  separate  parts,  a  part  r  (t)  significant  only  at  low 
t  and  a  part  r2(t)  significant  only  at  high  t.  Accordingly,  from 
(15) , a (T)  is  also  a  sum  of  two  parts, 


r:,f 


a  (T)  =  (T)  +  a2  (T)  , 


(16) 


each  defined  by  analogy  with  (15) . 

Representation  of  a, (T). 

For  convenience  we  define  a  Step  Function,  H(x),  of 
x  as  follows: 


and 


x  <  0 ; 
x  >  0 ; 


H  (x)  =  0 
H  (x)  =  1 


(17a) 

(17b) 


We  assume  that  the  mathematical  form  of  r  (t)  for  all  t  is  the 
same  as  that  of  the  ID  y(t)  at  low  t;  that  is. 


r,  (t,e,4>,x,y)  =  t  /t  -  (t  /t-1)  H(t  /t-1) 


(18) 


where  t  =  t  (9,<};,x,y)  in  general  for  a  particular  gas-solid 
system.  Thus  r.  is  significant  only  at  low  t,  as  required.  We 
consider  three  Assumptions  as  to  the  form  of  r  . 

Assumption  1:  t  is  independent  of  0,<f>,x  and  v  and  may  be  written 


t  =  T 
c  o 


(19) 


This  assumption  combined  with  (18)  and  (15)  gives 


ri(t)  =  "  (To/t  "  1]  H^T0/t"1) 


(20) 


and 


aj  (T) 


1  -  «xp(-  t2)  ~  2I*  exp  (-  |) 


(21) 


Assumption  2 :  t  depends  only  on  9 ,  and  this  dependence  is  such 
that  f l  depends  only  on  the  "normal  component",  kt  cos2 6,  of  the 
incident  energy: 


^4 


V; 
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tc  =  To  sec2e  (22) 

rT(t)  =  (To/t)sec2e  -  [ (T  /t) sec2 6-1] K[ (T  /t)  sec2 e-I]  (23) 

and 

T  T  T 

ax(T)  =  1  -  exp  (-  ;p-)  +  ^  eif  ~)  (24) 

where  ei(x)  is  an  exponential  integral15  of  x: 


ei(x)  =  J  exp(-U  d £/£  •  (25) 

x 

Assumption  3;  Each  "compohent"  of  the  incident  energy,  the 
"normal  component",  ktcos26,  and  the  "tangential  component", 
ktsin2e,  is  accommodated  according  to  a  formula  of  type  (18), 
and  each  "component"  is  associated  with  the  same  critical  initial 
effective  temperature,  t  ,  expressed  as  in  (19) .  This  assumption, 

v*» 

which  is  the  same  as  that  inherent  in  the  lattice  theory,  implies 
that 


r  ft)  =  2T  /t- (T  /t  -  cos2 6 ) H (T  /t-cos2 0 ) - (T  /t-sin2 0 ) H (T  /t-sin2 0 ) 
x  o  o  o  o  o 

(26) 


and  leads  simply  to 

oj  (T) 


(27a) 


We  note  that  the  relation  (27)  follows  also  by  combining  (18) 
with  (13;,  putting  a  =  a  ,  y  =  r.  and  t  =  T  .  For  low  T, 

A  •  v»  U 

L.T .  gives  a  similar  result: 

t 
T 


(low  T) 


a  (T) 


1  -  exp 


(27b) 
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None  of  these  three  Assumptions  gives  results  very 
different  from  either  of  the  others  in  the  ranges  of  interest 
(T  appears  as  a  dispensible  parameter) .  We  choose  (27)  ,  that 
is  Assumption  3,  to  represent  ot^T)  on  the  basis  of  the 
simplicity  of  the  result  and  the  fact  that  it  may  be  considered 
a  first  approximation  to  either  of  the  other  results,  and  in 
fact,  has  a  value  which  lies  between  them. 

Representation  of  a2 (T) 

We  let  be  a  "characteristic  frequency"  associated  with 
the  solid  surface,  and  an  "interaction  time"  of  the  gas- 
solid  interaction9-11.  We  assume  that  the  hard  spheres  limit 
is  approached  asymptotically  as  0,  and  a  first  guess  at 

the  form  of  this  approach  may  be 

r2  (t,e,$,x,y)  -» 

(large  t)  -  =  1  -  exp  (-  7 — *  (28a) 

r2  (»,e,<fr,x,y)  “Vi 

However,  we  find  that  the  experimental  data  are,  in  general, 

approximated  better  by  the  form 

r2  (t,  8  ,  <f>  ,x,y)  .. 

(large  t)  - -  =  tanh  (  —  -  .)  (28b) 

r2  (w,e  ,<j>  ,x,y)  ^o1  i 

and  we  use  the  form  (28b)  throughout  this  paper. 

We  assume  that,  at  least  for  small  values  of  y,  is 

related  approximately  to  the  normal  component,  UQ  cos 9,  of  the 

incident,  velocity  of  the  gas  atom  normal  to  the  solid  surface 

and  to  the  Morse  parameter  a  of  the  gas  atom-solid  atom  potential 

bv 


if* 


10 


at,U  cos6  ;  constant  (29) 

l  o 

Also,  Uq  is  related  to  t  via  (10a)  to  obtain 

r  (t,9,<J>,x,y) 

(large  t)  -  =  tanh[c'  —  (^— )  cose]  (30) 

r2(«,e,*,x,y)  “o 

where  c*  is  a  constant  of  order  unity. 

L.T.  is  concerned  with  evaluation  of  r  (t,0,4>  ,9,0)  and 
for  ti  <  *v»  0.84  the  hard  spheres  limit  of  this  function  is10-'11 

r  (“,0,4>,0,0)  =  4p/(l+u)2  .  (31) 

The  form  of  r  («,8,$,x,y)  is  obtained  through  consider¬ 
ation  of  results  from  a  3D  theory  of  surface  scattering12  cited 
above,  an  approximate  result  of  which  is 


A  !  F  (“,  0 ,  <j>  ,x,y) dx  dy  =  — cose  .  (32) 

J  (1+u)2 

Oman16  proposed  a  similar  formula  in  which  the  3.6  on  the  R.H.S. 
is  replaced  by  4.  However,  we  consider  (32)  more  nearly  correct. 
Combining  (32)  with  (30)  gives 


(large  t)  ~ 


r  (t,6,$,x,y)dxdy 


3 - cose  tanh[c'  —  cosO] 

(1+W)2  “o  M 

(33) 


where  we  have  used  the  fact  that  r2  (»,e  ,<j>  ,x,y)  =  r  («*>,  8 , $ ,x,y)  . 

We  assume  further  than  an  approximate  form  of  a2  (T)  for 
large  T  may  be  obtained  from  (15)  and  (33)  by  replacing  cose 
and  t  in  the  integral  by  their  respective  mean  values  via 
(14)  . 
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and 


co s 6  =  2/3 
t  =  2T 


(34a) 


to  give 


<x2  (T)  - 


_  2 . 4y 


tanh 


a  f  kT^I 


to  VM 
0 


(35) 


(1+y) 2 

where  c  is  a  constant  of  order  unity.  In  order  that  a  '  (T)  +  a2(T) 
should  not  exceed  unity  at  smaller  T  we  multiply  a by  exp(-To/T) 
to  obtain 


a  (T) 

2 


~  a  (=°)  tanh 


*1  r 


am  rTyz  j  r  o\ 

Tr  '-M''  t  J 


(36) 


where 

a(«)  =  2.4y/(l+y)2  (37) 


and  X '  is  a  parameter  whose  value  depends  only  on  the  solid 
surface. 

The  arguments  above  may  have  to  be  modified  for  larger 
values  of  y  because  of  the  possibility  of  several  collisions  at 
each  encounter  between  gas  atom  and  solid  surface.  In  principle, 
the  nature  of  these  modifications  is  not  clear.  In  practice, 
we  find  better  overall  agreement  with  experimental  data, 
especially  those  relevant  to  systems  with  larger  values  of  y , 
when  the  argument  of  the  tanh  function  is  multiplied  by  (1+y)? 
and  X f  is  replaced  by  a  new  X  whose  value  again  depends  only  on 


the  solid  surface. 
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-  - 
-s-V 

h-, 

J-4 


\i 


Our  representation  of  a,  (T)  is 

a2 (T)  =  a (“)  tanh 


(MT)^ 

aT~T 


(38) 


where  the  hard  spheres  limit  a{®),  given  by  (37),  is  assumed 
valid  for  all  values  of  p;  this  assumption  is  supported  by 
qualitative  evidence  from  the  lattice  theory8"11. 

Our  formula  is  obtained  by  combining  (16)  with  (27)  and 

(38)  : 


a (T)  =  1  -  exp( — +  a (»)  tanh 


(MT) 

ST 


'/* 


T 


a 

X 


exp  (-  =£) 


(39) 


The  constant  c  of  (35)  may  be  related  to  the  X  of  (39)  by 
equating  the  arguments  of  the  two  tanh  functions.  We  assume 

that  the  characteristic  frequency  of  the  surface,  w_  ,  is  the 

0 

Debye  frequency  of  the  solid;  that  is 


u>o  =  k0/fi 

and  we  obtain,  for  small  p. 


c 


me 

17X 


2 


mole  gm 


(40) 


(41) 


PROCEDURES  AND  COMPARISON  WITH  EXPERIMENTAL  DATA 


We  use  two  different  types  of  experimental  data,  one 
giving  the  dependence  of  a  on  T  for  several  gas-solid  systems, 
the  other  consisting  of  a  values  for  several  gas-solid  systems, 
each  determined  at  a  single  T  value.  We  select  only  data  from 
experiments  in  which  we  think  considerable  care  has  been  taken 
to  control  the  state  of  the  surface,  and  to  our  knowledge  we 
include  all  available  data  of  this  nature.  Discussion  of 
available  ac  data  and  the  consequences  of  failing  to  excercise 
reproducible  control  in  the  course  of  experiment  is  given 
elsewhere1.  We  use  one  of  two  procedures,  depending  on  the 
nature  of  the  data,  to  check  our  formula  against  experiment; 
these  two  procedures  (A  and  B)  are  described  below  together 
with  a  third  (C)  which  is  useful  in  special  cases. 

Procedures 

Procedure  A  is  used  when  temperature-dependent  data  are  available 

for  at  least  one  gas.  If  data  are  available  for  more  than  or.e 

gas  on  the  same  solid,  it  is  convenient  to  begin  with  the  gas 

which  a_  any  given  temperature  is  suspected  to  adsorb  least 

on  the  given  surface,  that  is,  the  gas  associated  with  the  lowest 

value  of  T  (He  in  all  cases  we  considered  here) .  We  then  find 

the  "best  fit"  values  of  T  and  X,  i.e.,  the  values  for  which  our 

o 

Formula  (39)  fits  the  T-dependent  data  as  well  as  possible.  The 
resultant  value  of  Tq  is  appropriate  to  the  gas-solid  system  and 


1  .1 
d.  *x 


that  of  A  to  the  particular  surface  under  consideration.  This 

value  of  A  is  used  with  the  T-dependent  data  for  other  gases 

on  the  same  surface  to  obtain  the  best-fit  value  of  T  appropriate 

c 

to  each  of  the  other  gas-solid  systems. 

Procedure  B  is  used  when  ac  data  on  a  given  surface  are 
available  at  only  one  temperature.  It  may  be  properly  applied 
only  when  the  ac  on  a  given  surface  is  known  for  at  least  one 
gas  for  which  Tq  is  small;  otherwise  the  calculated  value  of  A 
is  liable  to  considerable  error.  For  any  particular  'rface  we 
use  the  He  ac  whenever  available  to  find  the  value  of  A  for  the 
surface  (setting  Tq  =  0  for  He) ,  and  then  use  this  value  of  A 
to  find  the  value  of  TQ  for  each  of  the  other  gases  on  this 
surface. 

Procedure  C  provides  a  quick  method  for  finding  a  value 
(approximate  in  general)  of  A.  It  is  applicable  when  data  are 
available  in  a  temperature  range  where  (T)  as  defined  by  (27) 
is  negligible  compared  to  a{T),  and  where  the  difference  between 
the  tanh  function  in  (39)  and  its  argument  is  negligible.  For 
an  accuracy  of  5%,  the  first  conditions  implies  Lhat 


20T 

T  >  - -  (42) 

a  (T) 

and  the  second  that 


a  (T)  <  0.4  a  (°°)  .  (43) 

When  both  of  these  conditions  prevail  (39)  may  be  approximated  by 


a(T) 


(MT)‘/ir  a/ A  . 


(44) 
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Thus  one  exact  value  of  in  this  temperature  range  is  sufficient 
to  obtain  a  good  approximation  to  the  value  of  X.  In  practice, 

T 

measured  ac  values  may  be  plotted  against  T  to  yield  a  straight 
line  through  the  origin  whose  gradient  determines  X .  This  value 
of  X  may  be  used  with  experimental  data  at  lower  values  of  T, 
where  oij  (T)  is  not  negligible  compared  to  a (T) ,  to  find  Tq;  this 
value  of  Tq  may  then  be  used  to  correct  the  previously  estimated 
value  of  X,  and  a  new  value  of  TQ  determined. 

Comparison  with  experimental  data 

Most  of  the  data  on  the  temperature  dependence  of  the 
acs  of  the  inert  gases  on  tungsten  are  due  to  Thomas6  and 
Silvernail7;  the  others  (at  higher  temperatures  for  the  He-W 
system)  are  due  to  Goldstein  and  Tho-Nhan17  and  to  Waehman18. 

The  data  of  Refs.  6  and  7  were  obtained  with  both  simple- filament 
and  potential-lead  tubes.  The  data  from  simple  filament  tubes 
(represented  by  open  circles  in  Figs-  2-4)  have  been  corrected 
for  end-losses1'19  while  other  data  (represented  by  diamonds) 
are  uncorrected  from  potential  lead  tubes.  The  data  of  Ref.  17 
(squares  in  Fig.  2)  are  from  a  potential-lead  tube  while  those 
of  Ref.  18  (Triangles  in  Fig.  2)  are  from  simple-filament  tubes; 
neither  set  is  corrected  for  end  losses.  B’ailure  to  apply  end- 
loss  corrections  to  data  from  simple  filament  tubes  results  in 
ac  values  which  are  slightly  too  large  (about  8%  for  the  data 
of  Ref.  lb) .  All  the  experimental  tubes  except  those  of  Ref.  18 
contained  an  evaporated  mischmetal  getter  on  the  tube  walls  to 
aid  in  maintaining  clean  filament  surfaces. 
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Values  of  gas  pressures  used  in  calculating  acs  from 
experimental  measurements  may  require  correction  for  thermal 
transpiration  if  pressures  are  determined  when  the  McLeod  gauge 
and  the  experimental  tube  are  at  different  temperatures?  such 
corrections  have  been  applied  to  the  data  of  Refs.  6,  7  and  18. 

In  our  opinion,  the  only  temperature  dependent  ac  data  available 
on  clean  surfaces  are  on  tungsten  and  the  most  reliable  are 
those  of  references  6  and  7  in  the  temperature  range  77  -  303°K. 
Procedures  used  by  Raines20  to  obtain  the  ac  of  He  on  Ni  (see 
Fig.  5)  suggest  to  us  that  the  resulting  surface  is  not 
characteristic  of  the  clean  metal,  although  it  may  be  of  ne  ly 
constant  composition  over  the  experimental  temperature  range. 

The  experimental  ac  values  on  Na  and  K  from  77°  to  298°K 
obtained  by  Petersen21  were  determined  on  surfaces  of  the  bulk 
metals  (rather  than  on  monolayers  of  these  metals  atop  a  foreign 
metal)  which  were  presumed  to  be  free  of  adsorbed  gases. 
Blickensderfer 2 2  working  in  the  same  laboratory  after  Petersen 
found  that  alkali  metals  interact  with  glass  to  form  hydrogen 
and  concluded  that  this  gas  was  very  likely  present  during  the 
ac  measurements  cited  in  Ref.  21.  Petersen21  assumed  that 
hydrogen  does  not  chemisorb  on  bulk  potassium  at  room  temperature, 
and  cited  the  work  of  Trapnell 2 4  in  partial  support  of  this 
assumption.  However-  the  reaction  vessel  used  by  Trapnell24 
was  of  Pyrex,  and  the  work  of  Blickensderfer  suggests  that  hydrogen 
was  very  likely  present  during  these  measurements  also.  Indeed, 
Trapnell23  noted  the  presence  of  a  foreign  gas,  which  he 
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identified  as  hydrogen,  in  the  process  of  studying  adsorption 
of  C2H2  on  potassium;  however,  he  attributed  its  origin  to  a 

reaction  of  K  with  C2H2.  Accordingly,  the  potassium  surface 
used  by  Trapnell  might  have  been  hydrogen  saturated  before  being 
deliberately  exposed  to  the  adsorbate  gas.  Because  of  these 
considerations,  it  is  our  opinion  that  the  states  of  the  surfaces 
in  Petersen's  experiments  are  uncertain.  Neither  the  formula 
nor  the  lattice  theory  can  be  made  to  agree  satisfactorily  with 
the  temperature  dependence  exhibited  by  these  data;  for  these 
reasons  we  have  used  only  room  temperature  data  from  Ref.  21 
and  Procedure  B  on  the  assumption  that  at  room  temperature  each 
surface  is  of  constant  composition  in  its  interaction  with  He, 

Ne  and  A. 

We  use  twenty-one  ac  values  of  inert  gases,  each  at  one 
temperature  (see  Table  I) ,  on  eight  metal  surfaces,  all  determined 
in  Thomas '  laboratory  except  the  value  for  the  ac  of  Ne  on  Fe 
due  to  Eggleton  and  Tompkins30.  Another  set  of  ac  value*.- 
cite  below,  due  to  Wachman31  consists  of  data  on  the  ac  of 
He  and  Ne  near  room  temperature  on  initially  clean  W  deliberately 
covered  with  known  adsorbates. 

The  theoretical  curves  obtained  on  applying  Proceuure  A 
to  the  temperature  dependent  data  ^or  He,  Ne,  A,  Kr  and  Xe  on  W 
[Refs.  6,  7,  17  and  18]  and  for  He  -  Ni  [Ref.  20]  are  shown  as 
continuous  curves  in  Figs.  2-5.  The  dash-dot  curve  in  Fig.  3 
which  agrees  better  with  the  Ne-W  data  than  does  the  continuous 
curve  is  obtained  when  a2(T)  is  represented  by  (36)  rather  than 
by  the  modified  form  (38) .  Agreement  with  the  He-W  data  is 
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the  same  regardless  of  which  of  the  two  forms  for  a2  (T)  is 
chosen  because  Procedure  A  begins  with  the  He-W  data.  Hence, 
for  W,  the  A'  of  (36)  is  related  to  that  from  (38)  by 

X'  =  2.30  X  .  (45) 

Our  reason  for  selecting,  in  general,  the  form  in  (38)  over 
that  in  (36)  are  stated  above.  For  the  case  of  Ne-W,  we  display 
al (T) ,  c*2 (T)  and  a (T)  separately  in  Fig.  6,  using  the  form  (38) 
for  <*2  (T)  . 

Dashed  curves  in  Figs.  2-5  represent  the  best  fit  curves 
obtained  from  L.T.  as  it  stands  in  Refs.  10  and  11,  where  the 
fitting  procedure  is  described.  Forms  of  the  response  function 
of  the  lattice  model  used  in  L.T.  calculations  for  A,  Kr  and  Xe 
on  W  are  approximate1 ?/ 1 1  while  those  used  for  He  on  W  and  Ni  and 
Ne  on  W  are  more  exact32,  in  all  cases  the  lattice  model  is  that 
described  in  Refs.  10  and  11  (all  the  spring  constants  associated 
with  a  given  metal  have  the  same  va-.  ;e)  . 

Tables  I  and  II  contain  lists  of  the  systems  treated 
by  Procedures  A  and  B  respectively.  Appropriate  values  of  a. 
and  T  are  included  in  Table  II.  The  sources  of  data  are  given 
in  the  Tai-es,  together  with  the  computed  values  of  X  and  c  for 
each  surface  and  the  values  of  u,  a{®)  and  T  for  each  gas- 
surface  system.  The  values  of  A v  and  T  1  for  He-W  and  Ne-W 
in  Table  I  are  those  values  appropriate  when  o;^  (T)  is  represented 
by  (36) .  Estimates  of  Tq  for  a  given  gas  by  Procedure  B  are 
obtained  where  a  values  on  the  same  surface  are  known  for  that 
gas  and  for  another  gas  having  Tq  =  0  (i.e.,  He  in  all  cases 
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we  have  treated) .  In  estimating  Tq  we  arbitrarily  set 

a 1  3  0,  implying  Tq  3  0/  if  a?  is  found  to  be  within  10%  of 

the  experimental  value  of  a .  We  have  no  data  for  the  system 

He~Fe,  and  we  assume  that  T  =  0  for  Ne-Fe. 

o 

L.T .  is  applied  only  to  these  systems  for  which  ac 
data  are  available  at  several  temperatures*  Accordingly 
the  values  of  the  dispensible  parameters  of  the  theory,  0,/a 
and  D,  and  the  resulting  value  of  t  derived  from  comparison 
of  these  data  with  L.T.  are  included  in  Table  I.  The  values  of 
0/a  and  D  in  Table  I  for  Be-W,  He-Ni  and  Ne-W  are  not  the  same 
as  those  cited  in  Refs.  8-11.  The  differences  in  the  values 
associated  with  He-Ni  arise  because  this  system  was  previously 
treated  by  a  perturbation  theory9  which  is  not  used  here;  the 
differences  in  the  values  associated  with  He-W  and  Ne-W  arise 
because  certain  errors  in  values  of  experimental  data  used  in 
Refs.  9-11  have  been  eliminated  ir.  the  present  case.  Also 
included  in  the  Tables  are  relevant  values  of  0  and  a.  The 
values  of  a  in  Table  II  and  those  in  Column  4  of  Table  I  are 
calculated  from  literature  values  of  the  a's  appropriate  to  the 
gas  atom-gas  atom  and  solid  atom-solid  atom  interaction  by 
means  of  a  combination  ruls9.  Values  of  0  from  the  sources  cited 
are  used  to  calculate  c  from  (41)  and  to  derive  the  values  of 
a  appropriate  to  L.T.  (column  II,  Table  I)  from  L.T.  values  of 
e/a. 


The  value  0  -  380°K  is  chosen  for  W  rather  than  the 
value  0  =  330°K  used  in  Refs.  9-11  because  the  former  gives 
better  overall  agreement  between  a  values  derived  from  L.T. 
and  the  calculated  values.  The  value  6  =  320°K  is  chosen 
for  Ni  from  Ref.  28  rather  than  the  value  &  =  37^°K  from 
Ref.  29  because  the  former  gives  the  same  value  of  a  for  He- 


Ni  as  that  for  He-W 
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DISCUSSION 


The  formula  (38) ,  represented  by  the  continuous  curves 
in  Figs.  2-5,  agrees  well  with  the  ac  data  used,  which  cover 
wide  ranges  in  a  and  T.  The  best  agreement  is  with  the  He 
data  on  W  and  on  Ni;  we  include  the  Ni  data  notwithstanding 
our  reservations  regarding  its  quality,  on  the  assumption 
that  surface  conditions  are  constant  throughout  the  temperature 


range. 

As  stated  above  the  representation  of  «2 (T)  for  Ne-W 
by  (36)  (dash-dot  curve  in  Fig.  3)  instead  of  (38)  results 
in  better  agreement  with  the  experimental  data.  It  is  evident 
from  Figs.  2-5  that  agreement  between  the  Formula  and  L.T. 
(dashed  curves  in  Figs.  2-5)  is  very  good,  the  best  agreement 
being  for  Ne-W  when  a2(T)  is  represented  by  (36).  In  all  cases 
except  those  of  He  the  curves  from  L.T.  fit  the  data  slightly 
better  than  those  from  tne  Formula.  An  important  difference 
between  the  Formula  and  L.T.  is  in  at®)  which  has  the  value 
2.4y/(l+y)2  in  the  formula  and  4y/(l+y)2  (for  small  y)  in  L.T. 

As  a  consequence  of  this  difference,  L.m.  in  its  preisent  form 
must,  at  large  values  of  T,  result  in  larger  values  of  a  than 
does  the  Formula.  The  divergence  between  them  begins  to  become 
evident  not  far  above  room  temperature  for  Re  on  W  and  Ni  (Figs. 


2  and  4) .  Differences  between  values  of  t  from  L.T.  and  the 
corresponding  values  of  T  from  the  Fcrmula  arise  in  part  from 
the  difference  in  the  a (®)  values.  Although  in  no  case  is  t 


very  different  from  T 


(Table  I)  it  is  found  that  t  <  T  for 

G  O 
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data  to  the  left  of  the  minimum  in  a(T)  (A,  Kr  and  Xe  on  W) 
and  t  >  T  for  other  data  (He  on  W  and  Ni,  and  Ne  on  V?)  .  When 

C  w 

for  the  system  Ne-W  a, (T)  is  represented  by  (36)  the  resulting 
value  of  T  '  is  very  close  to  the  corresponding  tc  (Table  I). 

From  the  values  of  T  of  Table  II,  we  note  that  in  no 

o 

case  does  the  value  of  a2  calculated  from  the  Formula  for  Ne 

differ  from  the  experimental  value  of  a  by  more  than  10%.  This 

result  lends  considerable  support  to  the  Formula  especially 

in  View  of  the  large  ranges  of  the  relevant  parameters  (e.g., 

mass  ratio  y)  through  which  such  agreement  is  obtained.  From 

Tables  I  and  II  we  note  that  the  general  trends  in  Tq  values 

for  A,  Kr  on  Xe  are  reasonable;  i.e.,  Tq  increases  with  increasing 

gas  molecular  weight  on  a  particular  surface.  However,  the 

values  of  Tq  for  Li  seem  very  low  and  Tq  =  900°X  seems  unreasonably 

large  for  Xe-Mo.  This  last  Tq  value,  however,  may  be  too  large 

because  of  the  very  large  effect  which  a  relatively  small 

experimental  error  in  the  a  value  may  have  on  it;  e.g.,  a  15% 

reduction  in  the  cited  experimental  value  of  the  ac  of  Xe-Mo 

would  reduce  T  to  460°K. 

o 

Values  of  a  for  He  in  Table  I  calculated  from  L.T. 
are  too  low  in  comparison  with  those  of  the  other  gases,  and  we 
consider  this  entire  set  of  a  values  unsatisfactory  (see  the 
discussion  in  Ref.  9) .  The  values  of  a  are  likely  to  be  in 
error  because  L.T.  gives  too  high  a  value  of  a(ro),  and  as  a 
consequence  too  large  a  slope  of  the  a (T)  curve  just  to  the 
right  of  The  effect  is  most  pronounced  for  those  systems 
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for  which  data  <ue  available  in  this  range,  viz,,  He-W  and 
He-Ni.  Since,  in  general,  large  a  values  lead  to  large  slopes, 
agreement  of  L.T.  with  the  He  data  is  forced  by  selecting  an 
unrealistically  small  value  of  a.  Realization  of  these  facts 
suggests  a  method  to  be  presented  in  a  later  publication  for 
improving  L.T. 

Values  of  the  parameter  X  and  the  parameter  c  calculated 
for  each  surface  from  X  using  (41)  appear  in  the  Tables.  We 
have  no  idea  about  what  are  realistic  values  for  X.  However, 
it  is  remarkable  that  in  nearly  every  case  the  value  of  c  is 
very  close  to  unity,  the  notable  exceptions  being  c  =  0.2 
for  Li  and  c  =  2.4  for  Pt.  If  the  proper  value  of  c  for  each 
surface  is  indeed  unity,  then  we  must  suppose  that  the  value 
a  =  0.038  for  He  on  Pt34  is  too  high.  However,  if  we  were  to 
postulate  on  the  same  bases  that  the  value  of  c  for  Li  should 
be  closer  to  unity  we  could  have  to  make  the  assumption  that  tl a 
experimental  value  of  the  ac  of  He-Li  is  too  low,  which  is 
difficult  to  justify. 

Several  authors  (cited  in  Ref.  35)  have  noted  that  the 
ratio  of  the  ac  of  Ne  to  that  of  He  on  the  same  surface  is 
approximately  equal  to  the  square  root  of  the  mass  ratio  of 
Ne/He ,  i . e , , 


‘Ne 


He 


2.2 


(46) 


We  believe  that  this  result  is  fortuitous  in  the  sense  that 


room  temperature,  where  most  of  the  relevant  data  were  taken, 
is  in  the  range  where  (44)  applies  approximately  to  both  He 
ana  Ne  on  many  surfaces.  Accordingly,  using  (44)  and  the 


approximation  aHg 


§Ne  9ives 


°Ne<T> 


( 
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(47) 


Fig.  (7)  contains  a  plot  of  aNe (T) /aH  (T)  for  a  tungsten 
surface,  using  the  Fcrmrla,  for  0  <  T  <  2500°K.  It  is  clear  from 
tiie  plot  that  (46)  may  be  said  to  hold  over  the  entire  range 
of  experimentally  accessible  temperatures  above  abour.  250°K. 

It  is  interesting  that  (46)  holds  not  only  on  clean  surfaces 
but  also  on  gas  covered  surfaces.  For  example,  according  to 
Ref.  31  the  experimentally  determined  ratio  aNe/aHe  at  room 
temperature  is  2.73  on  hydrogen  covered  tungsten,  2.53  on 
deuteruim  covered  tungsten  and  2.47  on  oxygen  on  tungsten.  We 
are  not  yet  able  to  deal  with  systems  in  which  the  surface  is 
gas  covered  because  we  do  not  know  how  to  estimate  a(°°)  for 
such  cases.  However,  on  the  basis  of  available  data  it  is 
reasonable  to  expact  that  a  plot  similar  to  that  in  Fig.  7  would 
apply  to  these  systems  also. 

Perhaps  the  most  reliable  prediction  that  can  be  made 
with  the  formula  at  this  time  is  about  the  behavior  of  a(T) 
for  He3  on  tungsten.  The  value  of  c  for  a  tungsten  surface  is 
known  and  the  value  Tq  and  hence  a  (T)  for  He3-W  is  likely  to 
be  small,  therefore  a(T)  can  be  estimated  with  some  confidence. 
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Since  Thomas6  has  expressed  his  intention  to  study  this 
system  experimentally,  it  appeared  worthwhile  to  include  a 
comparison  (Fig.  8)  of  the  a (T)  curve  for  He3-W  with 
that  for  He4-W,  using  the  same  value  of  Tq  =  C,2°K  for  both. 
The  a (T)  for  He3-W  is  not  evaluated  below  50°K  where  a  large 
error  in  Tq  would  affect  the  results  significantly. 
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by  combining  Columns  5  and  9  of  this  table. 

:ation  of  procedure  A  where  o  (T)  is  represented  by  (36). 
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TABLE  II 


(a)  Not  stated  in  "Ref .  30.  . 

(b)  From  values  of  M  and  m  in  Ref.  25. 

(c)  From  equation  (37)  cf  tnis  paper 

(d)  Evaluated  from  the  combination  rule  and  the  value  of 

a  for  Ile-He  in  Ref.  9  and  the  values  of  a  for  Ni-Ni  and 
W-W  in  Ref.  26. 

(e)  Guessed. 

(f)  Guessed  as  in  Ref.  9. 

(g)  Frora  the  combination  rule  in  Ref.  9  and  the  data  for  the 
gas-gas  and  the  solid- solid  systems  in  Refs.  26  and  27. 

(n)  From  application  of  Profcedure  B. 

(i)  Estimated  by  Procedure  B. 

(j)  From  (41)  of  this  paper. 
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